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f"S i The paraxial (parabolic) theory of a near forward scattering of a quantum charged particle by 

a static magnetic field is presented. From the paraxial solution to the Aharonov-Bohm scattering 
problem the transverse transfered momentum (the Lorentz force) is found. Multiple scattering is 
considered for two models: (i) Gaussian 5-correlated random magnetic field; (ii) a random array of 
the Aharonov-Bohm magnetic flux line. The paraxial gauge-invariant two-particle Green function 
averaged with respect to the random magnetic field is found by an exact evaluation of the Feynman 
integral. It is shown that in spite of the anomalous character of the forward scattering, the transport 
properties can be described by the Boltzmann equation. The Landau quantization in the gauge field 
[ of the Aharonov-Bohm lines is discussed. 
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I. INTRODUCTION 



| The paper addresses the problem of quantum transport of a charge in an inhomogeneous static random magnetic 
c*j . field. In recent years, this or related problems have been met in a number of contexts in physics of 2-dimcnsional 
systems. For instance, in the composite fermion model of the Fractional Quantum Hall Effect, a (fictitious) random 
magnetic field is the environment which controls dynamics of effective charge carriers Q. One meets the fluctuating 
gauge fields in some models of high temperature superconductors , where the gauge field is the tool to impose the 
£SJ | constrain of no double occupancy in the t — J model ||. Besides, stochastically inhomogeneous magnetic field can 
^ • be experimentally created by various ways. For example, the field is irregular near the surface of a superconductor 
in an external magnetic field if the Abrikosov flux lattice is disordered; depending on the experimental conditions, 
the magnetic field inhomogeneities may be weak and smooth, or the field may be concentrated in an irregular array 
of flux tubes. Various aspects of transport in the magnetic field of the Abrikosov vortices, the weak localization in 
particular, have been studied in Refs. p|-|l0||. In recent years, the random magnetic field problem has been an active 
• subject area Refs. ]TT|-[l9l|. 

The formulation of the problem is as follows. A particle with the electric charge e and the mass m moves on the 
x — y plane subject to a vector potential potential A(A X , A y ) generated by a magnetic field b(x, y) = (rot A) . Two 
random fields models are considered in the paper. In the first one, the magnetic field 6(r) is a random Gaussian 
£h variable with zero average, (b) — 0, specified by the correlator 

(&(r)6(r')> = (^J -p5(r-r'), (1.1) 

O ■ 

where $o = is the flux quantum; the strength of the random magnetic field is characterized via the length C 
• *h . the meaning of which is that the magnetic flux through the area C 2 is typically of order of $o- The random field is 
assumed to be weak in the sense that C much exceeds the wave length A = h/p, p being the particle momentum. In 
another model [2C[-|2^| which is motivated by fractional statistics theories, the gauge potential is created by a random 
array of the Aharonov-Bohm flux lines. A system of the Abrikosov vortices (e.g. in the gate of a MOSFET transistor 
|j| p[) may serve as an experimental realization of the Aharonov-Bohm array if the particle wave length much exceeds 
the vortex (magnetic) size. 

In the random magnetic field case, the traditional approach of the theory of disordered systems [E3| meets difficulties 
on the very first steps. Indeed, the simplest object that is the single particle Green's function G(l, 2) = (ip(ri)ip*(r2)), 
is not gauge invariant and the physical meaning of its averaging with respect to the vector potential A generated by the 
random field is not clear. One may define a gauge invariant combination G(l, 2) = (ip(ri)ip* (7*2)) exp[i^ J c ^ dl • A] 

where the path C12 connects the points t*i and 1*2. Albeit gauge independent, G(l, 2) essentially depends on the choice 
of the path C12. With the point T\ 2 connected by the straight line, the field averaged G(l, 2) has been found in Ref. 

int. 
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Another problem is the diverging scattering rate i. For small scattering angles <fi, the differential cross-section 
behaves like i so that the scattering total cross-section is infinite. In other words, the life time of a state with 

the definite momentum is zero. The conventional diagram technique pljfl where - is assumed to be small compared 
with the kinetic energy, becomes questionable. On the other hand, it is known |24f| that in gauge-invariant response 
functions the self energy enters in combination with the vertex corrections and the divergence cancels out. In Ref. 
[H , it has been attempted to introduce a physically sensible gauge- invariant "single-particle time" r as a parameter 
entering the Landau level broadening. 

The main purpose of this paper is to develop a scheme which allows one to study the most singular part of interaction 
with random magnetic field, that is the near forward scattering. 

To pinpoint the physics behind the theoretical difficulties, consider first propagation of a plane wave. It is common 
in wave mechanics to analyze propagation in terms of wave fronts, i.e. the surfaces (lines in 2D) of constant phase. 
The property of the wave front line is that the probability current is locally perpendicular to the line. In the magnetic 
field, the phase of the wave function is ill-defined because of the gauge freedom. Nevertheless, one can construct a 
gauge invariant quantity x defined on a line, which in a limited sense plays the role of the phase: Given the wave 
function ip(r) and the vector potential A, the phase for the points s on a line S is defined through its differential 
as 

Tfl 

d X = m J^ (1-2) 

where j — ^3? (?/>* (7V — |^4) tp) is the probability current density. Provided the line S does not have self- 
intersections, x is an unique function of s. 

If X is a constant, i.e. j-ds = and j _L ds, the local current and the normal to S are parallel, so that 5 is a 
wave front. If x is a slowly varying function, A (x( s ) — x(0)) gives the local distance from S to the wave front passing 
through the point s = 0. 





(a) (b) 
FIG. 1. The wave front of a wave on a plane, (a) The wave front 5* is a line at each point of which the probability current is 
directed perpendicular to the line. The gauge invariant phase x is a constant along a wave front. On a gener al li ne S, chosen as 
a straight line in (b), the phase a function of the coordinate along the line s. Its variation is found from Eq.( |l.2| ). The physical 
meaning of x( s ) is that the phase difference Ax(s) = (x( s ) — x(0)) multiplied by the wave length A shows the local distance 
AAx from the line S to the wave front. 

Consider now how the random magnetic field affects the wave front upon its propagation. Take a state ip, for which 
the line x = is a wave front corresponding to the propagation in the positive ^-direction (j x ^0). To satisfy the 
requirement dx — 0, the wave function is 



ip(x = 0, y) = exp 



eh . 

1 — / dy A y (x = 0,y) 



along the wave front x — (choice of lower limit of integration is not important). 

To find the profile of the wave front having advanced from x — to a finite x, one can apply the usual eikonal-type 
approximation where the field affects only the phase of the wave function through the factor exp[#- J dl ■ A), dl || x 
being along the direction of propagation: 



ip(x, y) = exp 



px 



dy A y (0,y) + - j dxA x (x,y) 
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Integrating Eq.(L2), one finds the phase x{V'i x ) as a function of y for fixed x. For the phase difference Ax(j/i, J/2j x ) 
x(V2l x ) ~ x{yi\ x )i one S e ts after simple calculations: 



$ 



(1.3) 



where $(2/1,2/2;^) is the magnetic flux through the area enclosed by the path (0,j/i) — ► (0,2/2) - * ( x iU2) — * (^jj/i) 
(0,yi) (see fig. g. 




x=0 x 

FIG. 2. Propagation of the wave front magnetic field. As discussed in the text, one constructs a state for which the straight 
line at x = is a wave front. For a finite x, the profile of the front passing via the point x,y\ is controlled by the phase 
difference Ax = x(2/2; x) — xiVi'i z )- The difference is proportional to the magnetic flux through through the shaded area. 



A non-local character of the interaction with magnetic filed is clearly seen from Eq.( |1.3| ): the phase difference is 
controlled by the flux rather than the magnetic field in the vicinity of the particle trajectories. The non-locality is 
obviously of the Aharonov-Bohm type. 



Averaging with respect of the random magnetic field Eq.(l.l), one gets the variation of the phase difference 



(Ax) (yi,y2; x ) = - 



1 M ■ 1 2/1 - 2/2 1 



C 2 



Most notable feature here is that (Ax) 2 grows with the separation \y\ — 7/2 1 (cf. Ref. pi[ ). For points separated by 
the distance Ay, the random phase difference is of order of 1 when the wave front advances to Ax ~ C 2 / Ay. One sees 
that an infinite plane wave, for which Ay — > 00, looses its coherence immediately, whatever small is the propagation 
distance Ax. (See Section ^ for a more formal derivation.) These qualitative arguments explain the actual meaning 
of the zero life time and show that it is not an artifact arising e.g. due to a violation of the gauge invariance. 

To handle the anomalously intensive forward scattering, one needs a method suitable for a nonperturbative analysis 
of the small-angle multiple scattering. For this, the paraxial (parabolic) approximation p5| to the Schrodinger 
equation is chosen in the paper. The paraxial theory is applicable when the particle moves mainly in the direction 
of an "axis" and the momentum transverse to the axis remains always small. The paraxial approximation to the 
wave equation is most popular in optics where it gives a convenient description of light beams propagating in optical 
systems, their diffraction, focusing etc |2(| . Taking scattering and diffraction broadening on equal footing, the paraxial 
approximation is more generally applicable then the eikonal one (27) . 

To make the paper self-contained , the derivation of the paraxial approximation is outlined in Sect.||. The case of 
magnetic field is considered in Sect. II A where a scheme for description of scattering by magnetic field is suggested. 
The scheme is in a sense gauge invariant, gauge freedom revea ling itself only in the overall phases. As a limiting case, 
one recovers the well-known eikonal approximation (see Sect. II B). 

To illustrate usage of the paraxial approximation, a simple problem of scattering of a charged particle by the 
Aharonov-Bohm magnetic flux line is considered in Sect. [II. This (or equivalent) problem is of interest in a broad 



variety of contexts extending from the cosmic string theory [|28| to superfluids 29-31 (the Iordanskii force) and 
superconductors ]32] | where the scattering of excitations by quantized vortex lines controls the vortex dynamics. 
Although the exact solution to the scattering problem has been known since the original paper of Aharonov and 
Bohm ]33( ] (see also review p4| ] and references therein), certain controversy in the analysis and the interpretation of 
the solution still remains. Different opinions exist in the literature about the existence of the transverse force exerted 
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by the Aharonov-Bohm line or a superfluid vortex. On the basis of the left-right symmetry in the Aharonov-Bohm 
differential cross-section, the authors of Refs. |}5| and Rcf. |3(| have come to the conclusion that the line docs not 
exert any Lorentz-like force (translated as the Iordanskii force in a superfluid). Other authors, |3l],^(],^2|,[57 38 predict 
a finite force. Due to its simplicity, the paraxial solution allows one to perform a detailed analysis and resolve the 
controversy. 



It is shown in In Sect. IV, that the paraxial scattering theory becomes manifestly gauge invariant if formulated in 
terms of by-linear in i\) and ip* object that is the density matrix p. The evolution of the density matrix is given by a 
gauge invariant two-pa rticle Green function. 

As discussed in Sect. [V A, the paraxial 2D stationary equation with inhomogeneous magnetic field can be written 
as a non-stationary ID Schrodinger equation for a particle in a time-dependent electric field. This mapping allows 
one to present stationary solutions to the 2D magnetic field problem as the Feynman path integral for the effective 
ID problem. 

In Sect.^, the paraxial theory is applied to the model of <5-correlated random magnetic field. It turns out to be 
possible to evaluate the Feynman path integral and by this to find a (paraxially) exact expression for the two-particle 
Green's function averaged with respect to the magnetic field fluctuations. It is shown that the density matrix evolution 
can be mapped to the Boltzmann kinetic equation. 

In Sect. VI, another model is considered where the random gauge filed is generated by a random array of Aharonov- 
Bohm fluxes. The flux lines are randomly distributed in the plane, the flux of a line, <I>, is distributed with the 
probability The Aharonov-Bohm array may create an effective magnetic field B if p($) is asymmetric, p($) ^ 

p(-$). 



In Sect. VI A, the Boltzmann equation for charge subject to a Gaussian random magnetic field or f ield o f AB-array 
is derived. With the help of the Boltzmann equation, the resistivity tensor is found. Finally, in Sect. VI B we discuss 
the density of states of the levels due to the quantization of motion in the field B. 



The results are summarized in Section VII 



II. PARAXIAL APPROXIMATION 

The paraxial approximation allows one to construct a family of solutions to the Schrodinger equation which are close 
to to the plane wave with a certain momentum p . The wave function ^ scfn a solution to the stationary Schrodinger 
equation, is presented as 



(2.1) 



where the envelope paraxial function ^(r) is supposed to be slowly varying at the distances of order of the wave 
length X = 



Po 



The Schrodinger equation reads 



(E(p) + U)V 



Sch 



EV 



Sch 



(2.2) 



E(p) 



-p 2 and U(r) being the kinetic and potential energy respectively. Given po, the family of solutions in 



Eq.( [2.l| ) corresponds to the eigen-energy E — E(po) and the velocity v = d ^^f 1 ■ Inserting Eq.( |2.l| ) into Eq.( |2.2| 
one gets equation for 

(e(p) - E{p Q ) + [/)* = , E(p) = e-k p °- r E(p)e^ p °- r , 



which is still exact. The operator E(p) = E(p Q + | V) acting on the slowly varying function 'J is approximated in 
the paraxial theory as 



E(p) a E(p Q ) + -v-V + (?Vi 



here Vj^ denotes the gradient in the direction perpendicular to v. 
The paraxial approximation to the Schrodinger equation reads 



ihv-'V 



2m 



-VI - J7 * = 



(2.3) 



4 



(Condition of applicability are discussed later). The main feature of the paraxial approximation Eq^^), is that it is 
of first order differential equation relative to t he coordinate in the direction of the propagation x = r • v/v. 

Introducing formally "time" r = x/v, Eq.( [2.3| ) takes the form, of the time dependent Schrodinger equation in a 
reduced space dimension: 



• t d T ( h \r2 
ik—V = -—VI 

8t \ 2m - 1 - 



U * 



(2.4) 



This formal analogy allows one to discuss the stationary solutions in terms of the wave m ovin g in the direction v, 
and call x the current coordinate of the wave. The Feynman path integral equivalent to Eq.(2.4) gives an alternative 
method of solving the equation. 

The probability current J is derived from the standard expression J = — 3f\ 
In the main approximation, the components parallel, Ju, and perpendicular, Jj_, relative to v, are 



: ^* Sch V^ S ch and the definition Eq.fl2.1D 



'II 



(2.5) 



These expressions are consisten t w ith the current conservation and Eq.(2.3) or Eq.(2.4): Indeed, the continuity 
equation which follows from Eq.(2.4), 



~d7 



div J ± = 



(2.6) 



The continuity in Eq.(2.6) means that the paraxial wave function 
dr-J^I 2 = 1 



is equivalent to div J = with J from Eq. 
can be normalized: 



fixing the total flux in the beam to v. 

The required solution to the paraxi al eq uation can be chosen by imposing a proper boundary condition to Eq.( [2.3| ) 
("initial" condition in the case of Eq.(2.4)) 

a. Paraxial approximation: conditions of applicability The paraxial theory is based on the approximation E(p) — 
E(po) ~ v5p\\ + i^(Sp±) 2 , 5p = p — po where the term 2^r(^P||) 2 is neglected. This is justifiable if the angle, 9 ~ 

between p and po is small, 9 <C 1. If the motion is free, 5p± ~ ^ where w is the width of the beam (defined by 
the boundary conditions). Paraxial approximation is therefore applicable if the beam is wide in the scale of the wave 
length A, 



5p± 
Po 



» A . 



(2.7) 



The scattering by the potential U changes the angle by 9 scat ~ U/E . The paraxial approximation requires the 
small angle scattering to be dominant, 9 ~ sca t *C 1, so that the theory is applicable only for fast particles: E 3> \U\, 
It is important however, that, as in the case of the eikonal approximation | p7j , the theory is applicable beyond the 
Born approximation: the phase shift Sep ~ j^j ~ 'B^f ' a being the thickness of the layer where U ^ 0, may be large 
even for fast particles. 

Unlike the eikonal approximation where only the phase variations are taken into account, the paraxial theory allows 
also for the change of the profile o f th e beam, i.e. \^{r)\ due to diffraction (or, equivalently, to broadening of the 
wave packet in the language of Eq.(2.4)). If w is a typical size of the transverse structure defined by either the initial 
condition or scattering, the "diffraction blurring of the image" happens when the beam travels the distance XdiS, 



XdiS ' 



w 

T 



(2.8) 



The diffraction length x^is is the typical distance for the paraxial approximation while region x <C Xd\s is described 
by the eikonal approximation. 

Applicability of the approximation at large distances requires further analysis. The paraxial relation 5p\\ = 
-(5p±) 2 is valid up to a small correction (Sp\\)% due to the neglected quadratic term: (<5p||)2 — (^P\\) 2 - ^ n the 



l 

2po ' 



main approximation, (<5p||)2 ~ P±/Po- Although small in comparison with 8p\\, the correction is important at long 
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enough distances; it can be ignored only if (8p\\)2x/% < 1. Thus, the paraxial approximation is reliable if the distance 
traveled by the beam is not not too large, 



x < 



-r I XdiS 



(2.9) 



If the main condition Eq.(2.7) of applicability of the paraxial approximation is met, this requirement is distances large 
compared with the typical diffraction length XdiS- 



A. Paraxial approximation: magnetic scattering 

In this section the 2D paraxial theory is applied to the case of an external magnetic field; for simplicity U = 0, 
generalization to U ^ is straightforward. 

The paraxial wave equation, the gauge covariant form of Eq.( [2.3| ) with U = 0, reads 

(i^d x + ^)i; = (2.10) 
where d XiV = i-^A xy A being the vector potential; the a;— axis is chosen in the direction of the propagation v. 



The current density is given by Eq.(2.5) if modified by the standard diamagnetic term |27J. 

It is convenient to consider first the situation when the field is present only in a finite region. Divide the space into 
three regions: x < Xi n , incoming (I); Xi n < x < x ut scattering (II); and outgoing region (III), x > x ou t (see Fig. ||). 
In regions I and II a magnetic field is absent. Present the wave function in I as 

*f J dr A 

*(r) = e R ' Vm(r) , x < x in , (2.11) 



and by this define ip- ln . The integral in Eq.( 2.11 ) does not depend on the path of the integration if the latter is in the 
magnetic free region I. The overall phase of ip- ln depends on the choice of Rj and the gauge of A. It is convenient to 
put Rj on the I-II interface, Rj = (x = Xi n , y = y*), with somehow chosen y*. 



III 



FIG. 3. Magnetic scattering geometry. Initially, the particle moves in the field free region 7; tpi n (x,y) defined by Eq.(2.11) 
has the meaning of the incoming wave function in the gauge where the vector potential is zero in I. The particle is scattered 
by the magnetic field in the region // and then moves freely in the region III. In the gauge where the vector potential is zero 
in III, the outgoing wave is described by ^wt(r). To solve the scattering problem means to relate the interface value of the 
outgoing wave, ip ou t(y) = ipout(x ut,y), to that of the incoming one, ipi n {y) = tpin{xi„,y). 



The new function ip- ln has the properties of the wave function in the gauge A — 0: It obeys the free equation 

dx 2m dy 



ft? 

ihv lT + ^7r^)^n = (2.12) 
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and the probability current is given by Eq.(2.5) (without any diamagnetic term). 

The wave function of the incoming beam is the input to the scattering problem: It is assumed that the problem of 
a free propagation in I is solved, and the incoming wave at the I-II interface, ipin(f)\x=x in = tpiniy), is known. The 
normalization condition 



Win(y)| 2 = 1 



makes the to tal (c onserving) flux in the cc-direction equal to the velocity v. With the above choice of Rj, the wave 
function Eq.( [2.11 ) at the boundary of the region I reads 



if J dyi A y {x in ,y 1 ) 

y(x,y)\ x=Xm = e y ip- m (y) . 

where the integration is performed along a piece of the straight line x = Xj 



(2.13) 



The wave at x > Xi n has to be found from the paraxial equation Eq.(2.1C), solved with Eq.(2.13) as the boundary 
condition. The solution in the both, scattering (II) and outgoing (III), regions may be generally written as, 



X if / dvi A y (x iv ,,yi) 

dy' G R (x 1 y;x ln ,y')e y * V'in(y') , x > 



(2.14) 



where the Green function, G R (r,r') solves 

(ihvd x + — d£)G R (r,r') =iKv5(r-r') , G R = for x < x' . 
\ 2m y I 

Similar to the above consideration of I (sec Eq.( [2.li )), one defines in III a new function, ^outj by 



l f / dr A 
*(r)=e R '" ipout{r) , x > x out ; 



(2.15) 



(2.16) 



again, the whole path of integration must be in the field free region III; choose the initial point of the integration 
pat h Rin — (x out , y*) (or any other point at the II — III interface). Analogously to ipin in I, Vw(»") obeys the free 
Eq.( |2.12 ) and is fully defined i n the whol e reg ion III by Vw(y) that is the boundary value at the II-III interface, 
iponujy) = tp(xout,y)- From Eq.(2.16) and (2.14), we see that 



Vw(y) = / dy' g(x out ,y;x in ,y')ip in (y') 



(2.17) 



where the Green function Q is 



-if J dy 2 A v (x,y 2 ) if J dyi A y (x' \yi) 

g R (x,y;x',y') = e ** G R (x, y; x' , y') e »• 



(2.18) 



Eq.( [2.17 ) combined with Eq.( 2.18| ) relates the outgoing wave amplitude ipout to the incoming wave ipi a , and thus gives 
a general solution to the magnetic scattering problem in the paraxial approximation. 
Known Vw, one finds ip ou t{x,y) in the outgoing region /// as 



4>out(x,y)= / dy'G (x-x out ,y-y')ipout(y') , x>x out 



(2.19) 



with 



G R (x,y)=9(x) 



1 



\j2niXx 



(2.20) 



being the free propagator 
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B. The eikonal approximation 



If the scattering region II is narrow enough compared with the diffraction length Eq.(2.8), one can neglect the 
transverse derivatives in Eq.(2.10). This limit corresponds to the well-known eikonal approximation 27 1. The wave 
function obeys the eikonal equation (U — 0) 



ihd x ip — 



Solution to Eq.(2.15) reads in the eikonal limit 

if J dx" Ax(x",y) 

G R (x,y;x',y') = e " 



5{y-y') 6{x-x') . 



(2.21) 



(2.22) 



Substituting Eq.(2.22) into Eqs.( [2.1g ) and (2.17), one gets the eikonal solution to the magnetic scattering problem, 

^out(y) = e < jfe*<»tyin(i/) (2-23) 

where &(y) is the total magnetic flux through the directed area (see Fig. ||) restricted by the path (xi n) y*) — > 
(x in ,y) — ► (x out ,y)) —> {x out ,y*y, choice of y* is arbitrary affecting only the overall phase of Vw @- 




' out 



FIG. 4. The eikonal solution to the magnetic scattering problem. The scattering amplitude (see Eq.(2.2J)), is controlled by 
the flux encircled by the contour shown in the picture. 



III. ILLUSTRATIONS: AHARONOV-BOHM SCATTERING 



To illustrate the usage of the paraxial theory, the scattering on the Aharonov-Bohm magnetic line is considered in 
this section as an example. Some results presented in this section has been published in a short communication Ref. 

The scattering problem is formulated as follows. The incident wave ipi„ = ip- m {y) comes from the negative side 
of the cc-axis and the charge sees a magnetic line (extending in the z-direction), which creates the magnetic field 
b(x,y) — &o5(x)S(y). The final goal is to characterize the outgoing wave at x > 0. 

As the magnetic field is finite only at the origin, the scattering region II in the terminology of Section II A can 
be shrunk to just the line y = 0, so that Xj n = — and x ou t = +0- Seeing that scattering region is narrow, 



the eikonal approximation discussed in Section II B is applicable, except, perhaps, in the immediate vicinity of the 
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singular magnetic line point y = 0. The flux function, $(y), in the eikonal expression Eq.( 2.23 ) is easily found to be 
$(y) = $0(y) (if is chosen at -co) . The outcoming wave reads 



^out(y) = exp (-iir$sign(y)\ ipi n (y) 



(3.1) 



where $ = ^- Same result can readily obtained |3lj by solving the first order differential equation Eq.( 2.21 ) 

with the vector potential in the gauge 



A x (x,y) = --<$> sign(y)5{x) , A y (x,y) = 



The freely propagating outgoing wave is found from Eq.( 2.19j ) with x out = and i/w from Eq. (|3.l[ ) Ref. |37) , 

ipout(x,y) = ipm(x, y)cos7r$ sin 7r$ , a; > , (3.2) 

where ip- lrL {x,y) is the incoming wave continued to the region x > 0, 



ipin(x, y) = / dy'G£(x, y - y')ip in (y') 



(3.3) 



i.e. the wave in the absence of the Aharonov-Bohm line, and 



V(x, y) = / dy'G«(x, y - y') y> Vin(y') 



(3.4) 



Eqs.(3.2 3.4[ ) give the solution of the Aharonov-Bohm problem in the paraxial approximation for an arbitrary incoming 
wave ip in (y)- 

Using Eq.(pT20|), one can check that 



dy|Vin(x,y)| 2 = 1 , / dy\V(x,y)\ 



dy ip* n (x,y)V(x,y) = 0; 



for any x > and ar bitra ry (normalized) V'in; these relations guarantee the current conservation |42| . 

The solution in Eq.(3.2) gives a convenient tool for studying interaction of beams (waves packets) with the Aharonov- 
Bohm flux line. Some examples are considered below. 



A. Plane incident wave 



Plane wave corresponds to ipi n (y) = 1- From Eqs.(3.3) and (2.20), one gets ip- m {x,y) = 1, and from Eq.(3.4) Ref. 



V(x,y) = K(y) , y 



'2Xx 



(3.5) 



where K(t) is a familiar function, 



K (t) = JL I dt' e lt ' 2 



(3.6) 



related to the Fresnel integrals. The outgoing wave Eq.(3.2) reads 
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ipout(%,y) = cos 7r$ + iK(y) sin 7r$ 



(3.7) 



The complex wave function Vwt^, y) is a function of the real parameter y only. This means, that ipout(x, y) for any 
point of the half-plane x > spans a line on the complex plane $tip ~ Sf^>. This line (Fig. ||) is the Cornu spiral 
well-known in the diffraction theory [jl3|. The property of the Cornu spiral is tha t (d t) 2 is pro port ional to (dt) 2 where 
dt is the distance between the points corresponding to t and t + dt. From Eq.( |3.5|) and Eq.( |3.6|) , the element d£ of 
the leng th along the spiral (the arc-length) and dy are related as {dt} 2 = ^ sin 2 (7r<l)(<ii/) 2 . Since the wave function 
Eq.( |3.7[ ) is real for y = 0, the arc-length counted from the spiral point with Qifj = is proportional to y 



4>=1/16 

■!> 1 l 

4>=l/2 
<t>=3/4 
<K=15/16 




FIG. 5. The paraxial wave-function Eq.(|3.7|) on the ip -complex plane for the fluxes jg&o, j&o, f $0, |$o, if^o- The 
points on the spirals are uniquely characterized by the spiral arc- length counted from the point where Im ip = and measured 
in units proportional to sin$. For the point with the coordinates x,y, the wave function ip(x,y) corresponds to the point on 
the Cornu spiral with the arc-length equal to y = y/V2Xx . In accordance with Ref. the exact Aharonov-Bohm solution 
can be also mapped to the Cornu spiral; for a general scattering angle the mapping condition is that arc-length « yj(r — x)/X 
in proper units. 



The singular Aharonov-Bohm line problem gives a rather tough test for the validity of the paraxial approximation. 
The evaluation of its accuracy has been made simple by the recent observation Ref. 44 that the exact solution can 
be presented in the form Eq.(3.7) with y — > Y with Y found from 



Y 



(r-x) (\ - 



(3.8) 



r and (p being the cylindrical coordinates; the terms not shown in Eq. (|3.8[ ) are of order O (-) [Q. (Unexpectedly, the 
Aharonov-Bohm wave function on the x — y plane can be mapped on the Cornu spiral not only in the forward direction 
but for any scattering angle (and r>^), see Ref. 44 for details.) The paraxial approximation gives correctly the 

leading term < / « y in the vicinity of the forward direction and short wave length A -C r. 



B. Finite size beam 



The scattering of a plane wave by the Aharonov-Bohm line is highly singular: Eq.(3/7) shows that in the forward 
direction, ip <C i.e. y <C 1, the wave function equals to cos7rl> and does not converge at large distances to the 

plane incident wave as assumed in the standard scattering theory. Known from the exact solution |45|, this anomaly 
is the reason why the text-book scattering theory fails: The incoming plane wave and scattered wave cannot be 
separated ]45| ] and the scattering amplitude cannot be introduced as the object carrying the complete information 
about scattering |37|| . 

The singular behaviour is due to the combination of two factors: (i) the infinitely long range of the interaction with 
the line; (ii) the infinite extension of the plane wave. Since any physical state has a finite transverse extension W , the 
behaviour of the potential beyond the width W is irrelevant, and the singularities are expected to be regularized. 
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To show this, suppose that the incoming wave is beam-like with the profile 



where W is the beam width; it is assumed that W 3> X so that the paraxial approximation is applicable. The beam 
has a small but finite angular width <po = X/W. One can easily see from Eq.(|3~2]) that at large distances x 3> W 2 /X, 
the outgoing wave behaves like a spherical wave i.e. \ip ou t{x,y)\ 2 « P(ip)/x where ip = y/x is the scattering angle 
and Pi}f)i 



2X (<p sin 7r$ — ipo cos 7r<E>) 2 

+ vl? 



= — —2 1 . 2^2 ' ( 3J ) 



has the meaning of the angular distribution of the intensity in the outgoing wave. As expected, the distribution shown 
in Fig. is perfectly smooth. 



FIG. 6. Scat tering of a beam by the Aharonov-Bohm line. The solid line shows the angular distribution of the intensity 
P(ip) in Eq.fejj) for t he flux $ = +1/4 and the beam angular width </?o = 0.1. The dashed, line which is the Aharonov-Bohm 



cross-section Eq.(3.1C) fits well the distribution at large enough angles S> po where the intensity can be attributed to scattering. 



The central, \p\ < po, peak has a nontrivial structure, asymmetric relative to ip — > — p. The small angle left-right asymmetry, 
which is absent in the scattering cross-section, is responsible for the effective Lorentz force exerted by the Aharonov-Bohm line. 



For the angles larger than the beam angular width, i.e. \ip\ ^> ipo, one gets from Eq.( |3.9| ) that 



„, , 2 A sin 2 7r$ 
P(<P) « — , 

7T p z 

recovering the small angle asymptotics of the Aharonov-Bohm scattering cross-section [[33| 

d<r\ X sin 2 7r<l 



dp J 



AB 



2tt 



sin * 



(3.10) 



(3.11) 



It would be rather trivial if the angular broadening of the incident wave led just to a regularization of the forward 
scattering singularity. Most important is that a qualitati vely new feature becomes seen: Unlike the Aharonov- 
Bohm cross-section Eq.(3.11) the angular distribution in Eq.(|3.S|) is left-right asymmetric. The antisymmetric part is 



concentrated in the forward direction \p\ < po, i-e. within the angular width of the incoming wave i.e. in the region 
where the scattered and incident waves cannot be separated p6[ . The asymmetry means that the beam is deflected 
by the Aharonov-Bohm line as a whole |37],[38| ■ The order of magnitude of the deflection is the initial angular width 
<Po ~ X/W . One may say that the Aharonov-Bohm line not only scatters the incident wave but also modifies the 
unscattered part of the wave. Later, it will be shown that in a system of many Aharonov-Bohm lines, the deflections 
by individual lines coherently add together, and the lines act as an effective magnetic field 
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The deflection of the beam Atp can be presented via the momentum Apj_ = lA(/j transfered to the charge in the 
direction perpendicular to its initial velocity. A calculation details of which are collected in Appendix ^ gives the 
following result Ref. |37) : 

Ap ± = ft|V>in(0)| 2 sin27ri- (3.12) 

here the momentum transfer is expressed via the value of the normalized incoming wave at the position of the line. 
By comparison with the exact theory, the validity of this paraxial result has been confirmed by Berry ]3q| . 



IV. DENSITY MATRIX 



T he th eory of magnetic scattering presented in Section [I A is gauge invariant only in a limited sense. Although 
Eq.( |2.17 ) holds in arbitrary gauge, each of the objects there, i/'im Vwti and Q is gauge dependent, although only 
through the overall phase. For example, under A — > A + Vx simultaneously with \l/ — > e l "^ x , the incoming wave ip- lTL 
Eq.(2.11) is modified as — > e % i x ^ Rl ^ipi n ; we see also that the overall phase depends on the arbitrarily chosen i?/. 
Of course, the observables are independent from the global phase and the above gauge dependence does not create 
any problem. 

A truly gauge invariant theory can be formulated in terms of the by-linear in W and "density matrix" p(y\ ,y 2 ',x) 
defined as 



-if J dy' A y (x,y') 

p(y u y 2 ;x) = e " 2 V(x, yi )V*(x,y 2 ) 



(4.1) 



The "density matrix" p carries the full quantum information needed to find observables and is gauge invariant. In 
the field free regions, the density matrix is built from ipin and ip? n or ?/w and if)* ut ; these combinations depend on 
neither gauge nor Riji- 

The current density the x- and y- directions are (cf. Eq.([2~5|)) 



J x (x,y) =vp{y,y;x) , J y (x,y) = — - — ) p(y u y 2 ; x 



d 



d 



Seeing that the evolution of the wave function ^ in Eq.(4.1) is given by the propagator G Eq.( 2.15 ), the density 
matrix evolves from x' to x (x > x') as 



<yi=y2=y 



p(yi,y2-,x) 



dy[dy 2 g R (y 1 ,x;y' 1 ,x l ) piy'^y'^x^Q (y' 2 , x' ; y 2 , x 



(4.2) 



where Q R is defined by (2.18) and the advanced Green function Q A 

G A (rx,r 2 ) = (G R (r 2 , ri ))* . 
Introducing the two-particle Green function 

fC(yi,y2;x\y' 1 ,y 2 ;,x') = g R (y 1 ,x\y / 1 ,x')g A (y 2 ,x'\y 2 ,x) , 



Eq.(4.2) can be written as 



(4.3) 



p(yi,y2]x)= / dy' 1 dy' 2 IC(y 1 ,y 2 ;x\y' 1 ,y 2 ;,x')p(y' 1 ,y 2 - 1 x') 



(4.4) 



As before, the incoming wave enters the scattering problem as the b oun dary condition at x = Xi n : p(y\, y 2 ] Xi n ) 
Pin{yi, y 2 )- Further propagation of the incoming beam is given by Eq.(4.4) with x' = Xi n . 
For future references, we note the following property of the Green function: 
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-oo 
oo 



dyG R (x,y;x / ,y 1 ) G A (x\ y 2 ; x, y) = 9{x - x')5(yi ~ y 2 ) 
dyG R (x,y 1 ;x',y) G A {x ,y;x,y 2 ) = 0(x - x')8(y 1 - y 2 ) 



this relations expresses th e cu rrent conservation |47[ . 

In particular, from Eq.(4.5), one gets the conservation of the total current in the beam: 



dy p(y, y; x) = / dy p(y ,y \ x = 0) 



(4.5) 
(4.6) 



A. Path integral representation 



Similar to Eq^^), the stationary paraxial equation in 2D may be mapped to a time dependent ID problem: 



(4.7) 



where "time" r = ^, a — A y , and ip = —^A x . In the effective ID problem, the particle moves in the "electric field", 



F 



1 

-—a — vtp , 



(4.8) 



defined by the "vector potential" , a, and the "scalar potential" , ip. The gauge transformation, A — > A+ Vx, translates 
to a — > a + Vx, y — > <p — ^x, so that the effective electric field Eq.( |4.§| ) is indeed gauge invariant. Locally, the electric 
field is related to the magnetic field of the original problem b(x, y) as F — ~vb. 

The mapping to the effective ID problem allows one to use a convenient path integral representation for the 
paraxial propagators. Obviously, G R is just the retarded Green function for nonstationary equation Eq.(4.7) and in 
the Feynman path integral representation 

v 

G R (x,y;x',y') = J V[y{x)]e iS ^ 
y' 

where the action of the effective ID problem S = i J dT(my 2 /2 + eay/c — ecp) translates as 

X 



dr • A 



(4.9) 



y=y(x) 



where y x = . The action corresponding to Q R differs from Eq.(4.9) in the path of integration which should be 



extended in an obvious way to include the additional exponential fac tors in the definition of Q R Eq.(2.1 
With the help of Eq.(p~9|), the two particle Green function Eq^O) can be presented as 

K(yi,y 2 ;x\y[,y 2 ;,x') = J V[ yi {x)]V[y 2 {x)]e iS ^^ , 

where 



(4.10) 



S = So + ^$(fe/i],[j/ 2 ]) 



Sq being the free motion contribution, 
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So[yi(x),y 2 (x)] = — \ dx (yj x - y 2 2x ) , 



(4.11) 



and $ is the flux 

$([y 1 ],[y 2 ])= j> dr. A (4.12) 
threading the (oriented) area bounded by the paths y\ (x) and j/2 (x) and the vertical lines at x and x' (see FigJ^) . 




FIG. 7. The int egrat ion in Eq.( 4.10) is performed with respect to the paths yi(x) and 2/2(2;). Due to the phase factors 



introduced in Eqs.(2.11), and (2.16), the action contains the vector potential integrated along the closed loop shown here 



The path integral representation is used below to perform averaging with respect to the gauge field. 



V. GAUSSIAN RANDOM MAGNETIC FIELD 



This Section con cern s the averaging the two-particle Green's f unct ion Eq.(4.3) with respect to the Gaussian random 



mag netic field Eq.(l.l). In the path integral representation Eq.(4.10), the random field enters via the flux $([yi], \yi\) 
Eq.( 4.12 ). For the Gaussian field Eq.(l.l), the averaged value 



exp ( 2ni 



exp 



-A no {[yi 



2C 2 



where A n 



X 

>([yi - J/a]) = J dx I J/i - J/2 1 



J/2 J 



is the non-oriented area bounded by the paths y\{x) and 2/2(2:) and lines x = x' and x = x (see FigJ^). Further 

calculations are rather simple thanks to the fact that ^e l27rT o^ is a functional only of y-y — y 2 - This important 

simplification is a property of the models with a (5-correlated magnetic field. 
In the variables 

1 . 

y = Vi - J/2 , Y = - (yi + y 2 ) , 



14 



the kinetic energy contribution Sq Eq.( l.ll ) reads after integration by parts, 



So[yi(x),y 2 (x)] = -^VxY\%, ~ ^ / 



Since Y(x) enters only the 1S0, the integration e with respect to Y gives S(y xx ). This means, that the integration 

,2 

with respect to y(x) is limited to the path with y xx = %3L = that is the straight line connecting the initial and final 
points. After this, the integral is easily calculated. 



Finally, /C av that is the paraxial two-particle Green function Eq.(4.3) averaged with respect to the fluctuation 
magnetic field, reads 



K. w ,{y 1 ,y 2 ;x\y' 1 ,y 2 ;,x') = JC (yi, y 2 ; x\y[, y 2 ; , x') exp 



A 



2C 



(5.1) 



here /Co is the free two-particle propagator, and A is the (non-oriented) area formed by the straight line trajectories 
Figf 




FIG. 8. The straight line trajectories contributing to the Feynman integral. _4 no in Eq.(5.1) is the geometrical area inside 
the closed loop. If (j/i — 2/2X2/1 — y'2) < 0, the trajectories cross each other, and the area is built of two triangles. 



A = — x') x 



\Vi +Vi-V2 

(yi-y2) 2 + {y' 1 -y' 2 ) 2 
\yi-y-2\+\yi-y2\ 



y'2 1 



if (yi -mWi -y'2) > 
if (vi-mWi-y^) <o 



In the variables y, Y and (x — x') — » x, an expanded version of Eq.(5.1) reads 

1 



K av (y,y;y',r';x) 



2ttxX 



exp 



xX 



(Y-Y')(y-y') 



1 x 
8£2 



( 



V 



(y + y') 



IY2 



y 



y 



(5.2) 



The two-particle Green function in Eq.( p.2|) allows one to find averaged over the random field "evolution" of the 
density matrix and thus describes correlations in the gauge invariant observables like the density or t he c urrent. For 
instance, it describes transmission through a slab with a random magnetic field. In what follows, Eq.( ^3.2j ) is applied 
to some simple cases: (i) focusing of a coherent wave; (ii) the scattering of the partially coherent spatially uniform 
incoming wave. 
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A. Coherent propagation: Focusing 



Let the incident wave ipi n (v) be a converging Gaussian beam i.e. 



My) = , , 1/4 /r. e ' V( ^ + ^ ) . ( 5 - 3 ) 



(7r)!/4^ 

where / is the distance to the focal point and w <C / is the width of the beam. It is well-known that the "spherical" 



wave like that in Eq.(5.3) will converge at the focal point x = / producing a diffraction limited spot with the waist 
~ X/6 where 6 = w/ f <C 1 is the angular size of the beam as seen from the focal point. 

The distributi on of the averaged intensity for the beam propagating in a random magnetic field can be found with 
the help of Eq.(pT^). Taking for simplicity only the points on the beam axis y — 0, the intensity I{x) = p(0, 0;x), 
reads 

oo 



At the focal point x = /, 



I(x) = ^— / dY dye 
\ irw Xx ' 



I (f) = JJ I <>!>< 



In the limiting cases, 




x w , fw-^iC 2 



For the conditions when the magnetic field is not important, the upper line gives the usual diffraction limited value of 
the intensity; the width of the spot at the focal plane (line) is of order of ~ 1/I(f) ~ X/9, where 8 = w/f. The larger 
the aperture w, the larger the intensity in the focus and the smaller the size of the spot. However, in the presence of 
the random magnetic field, the intensity saturates when the aperture w ~ C 2 / f and 9 ~ (£/ f) 2 . 

This behaviour is very different from that when the scattering is due to a random scalar potential with a short 
correlation length. In this case, the relevant parameter characterizing the disorder is the focal length / in units of the 
mean free path I rather then the size of the aperture: On the background created by incoherent scattering, one would 
see a spot with the disorder insensitive profile and the integral intensity oc e^^ 1 (the exponential factor is probability 
that the wave does not experience any scattering). 

B. Incoherent wave: the Boltzmann equation 

Consider the initial density matrix of the form p^ n (2/1,2/2) — Po (j/i ~ V2 ) , V — y\—yi, which corresponds to a par tially 



coherent spatially homogeneous state. At x > 0, the density matrix p(yi, yi] x) = p(y; x) is found from Eq.(4.4) with 



the two-particle propagator from Eq.(5.2). The Fourier transform, 

00 

P{y,x)= I d(pe l(p % n v (x) , (5.5) 



defines n v which has the meaning of the distribution function with respect of the transverse momentum, q = ipj, 
and ip is the angle between the velocity of the particle and the x-axis. In the paraxial situation, the distribution is 



concentrated at s mall a ngle s and the integration in Eq.(5.5) can be taken in infinite limits 



From, Eqs.(4.4) and ( p.2| ) we easily get the density matrix of the wave having traveled the distance x: 

* I y 
2£ 2 



p(y;x) = p (y)exp ( — —-r- ) . (5.6) 
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As required by the current conservation, p(0; x) does not depend on x whereas the non-diagonal elements of the 
density matrix y ^ decay to zero, the faster, the more "distance to the diagonal" \y\. In other words, random 
magnetic field is very effective in destroying a long range coherence, the longer the coherence, the faster it decays. 



Considering the limiting case of plane infinite incident wave n v (x = 0) = 5(tp) i.e. Pi n (y) = 1, Eq.(p.6[) gives 



n v (x) 



7T If* 



A 2 



2C 2 



(5.7) 



The evolution is nonperturbative in the sense that the plane wave looses its shape immediately at any x ^ trans- 
forming into the Lorentz distribution with the width of the distribution proportional to x and the strength of the 
field |48[| . In particular, it means that the plane wave is not a good basis for the perturbation theory. 

More insight can be gained if the evolution of the density matrix is mapped to a Boltzmann-type kinetic equation. 
For this, note that the density matrix in Eq^ESJ)) satisfies the equation, 



(5.8) 



Written for the distribution function n v introduced in Eq.(5.5), the equation acquires the familiar Boltzmann form 



dx 



In,, 







(5.9) 



where I is the collision integral i.e. operator / in the ^-representation. 
One may present the collision integral in the standard form 



In v = / d(j) w{4>)(n tp - n v+t p) 



(5.10) 



where w(<f>) is the scattering rate for the process tp 



From the condition that the operators in Eq.(pT9|) 



and Eq.( |5.10|) have same eigenvalues corresponding to the common eigenfunctions, 
requirement that Ref. E9| 



w must satisfy the 



d<fiw(<fi)(l — e" 



(5-11) 



From here, 



ttt <p z r 4,L z 



(5.12) 



Usually, one can split the collision integral into the in- and out-scattering pieces. In the case of random magnetic 
field, the scattering-out rate is ill-defined as J d4>w((p) diverges at small angles, and the split hardly makes sense. 
On the other hand, the collision integral, as an operator acting on the distribution function, is well defined and the 
transport is not singular. 

Treating the random magnetic field in the Born approximation, Aronov et al. fl3|| found the the scattering rate 
W(<f>) to be 



W{cj>) 



1 



27TT 



■cot" 



(5.13) 



Because of divergence at small angles, one may doubt the validity of th e Born approximation. However, the small 
<j) asymptotics of W(tf>) agrees with Eq.( 5.12 ). This means that Eq.( 5.13 ) is actually v alid for arbitrary </> if used for 
constructing the collision integral. The other way around, the collis ion integral Eq.(5.10) is expected to cor rectl y 
describe scattering with arbitrary scattering angles if W(<f>) Eq.( 5.13 ) is used instead of paraxial w{<j>) in Eq. (|5.12 ). 
This allows one to generalize the paraxial kinetic equation, including large angle scattering. 

The kinetic equation for the distribution function n v reads 



v • Vr 



In v = 



< If < 27T 



(5.14) 



where the collision integral 
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2tt 

f 1 f d(f> 2 4> 
ln v = — J 2^ cot 2 ^ ~ rbip+ 't'> ■ 
o 

As before, the spilt of the collision integral in the in- and out- scattering parts leads to divergences and, therefore, 
has a very limited sense. At the same time, the collision operator is well-defined: If n v is presented as the sum, 

m 

over the eigenfunctions of the collision operator e lmv , m — 0, ±1, . . ., the collision operator acts as 

oo 

In v = - V(2|m|-l)n m e"^ , 

The parameter tq has the meaning of the relaxation time for the first harmonics m = ±1 i.e. the transport relaxation 
time. 

One concludes that transport of a charge in a random magnetic field can be described by the Boltzmann equation, 
and is not anomalous in spite of the fact that the total scattering rate is infinite. 

VI. RANDOM ARRAY OF AHARONOV-BOHM LINES 

This section deals with the model of a random gauge field where the gauge field is created by an array of Aharonov- 
Bohm lines. It is assumed that the lines in the array take random space positions and the flux of a line $ may be 
random. The model is specified by the averaged density of the lines g? ab and the probability distribution for the 



magnetic flux <f> in a line. Our primarily goal is to average the paraxial two-particle Green function Eq.(4.3) over the 
distribution of the lines. 

The calculations turn out to be very similar to those in Section M so that we only outline them. Repeating the 
arguments from Section ^ one comes to an expression similar to EqT ([5.l|) : 

r < , , i\ it f i ' >\l (<■> ■^ {t) (vuV2]x\y' 1 ,y2;,x') 

^a.v{yi,y2;x\y 1 ,y 2 ;,x ) = JC {yi,y2;x\y 1 ,y 2 ;,x ) < exp 2m- 



$0 

where $0 (yi, y 2 ; x\y[, y 2 ', , x') is the flux through the oriented area bounded by the straight (directed) lines connecting 
the initial and finite points (see Fig.||). Given the configuration of the Aharonov-Bohm array, the flux through the 
area is the sum over the lines piercing the area. The k's line with the flux <E>fc contributes to the total flux as <Jk&k 
where <7fc = +1 or — 1 depending on the orientation, positive or negative, of the area the line is situated in. Let 
N + (iV_) be the (random) number of lines in the area with positive (negative) orientation; The variables <!>fc's are 
independent in the model, and the averaging exp (27ri-|^-) over the configurations with fixed N± is simple: 

exp % — ) a k $k \ } = { exp 2tti — ) ( exp -2tti — 



, , yj K i K ■ , * v-ii-jj i Mil v if \ v^vjj i v 

$ o fr[ )/ \ \ ®oJ I \ \ *o 

where < exp(±^p) > implies averaging with the distribution function p($). 

The random numbers N± obey the Poisson distribution = e N N N /N\ with N either (N + + N-) — c? AB ^4 no or 
N + — JV_ = d AB ^4 , A no and A being non-oriented and oriented area, respectively. 
Finally, 

exp ( \ = exp (_jA*o + ^BA Q 



where A no is the non-oriented area Eq.(5.2), Ao is the oriented area 



A = --{x- x')(yi + y[ -y 2 - y 2 ) 
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and the effective magnetic field 



$ / $ 
B = d AB — ( sin 2ir— 

AB 27T \ $ 



(6.2) 



in Eqs.(6.1) and ( |6.2| ), the averaging is performed with the distribution function of the flux in the line £>($)• 
Collecting the results together, the two-particle Green function reads 



K av (y,Y;y',Y';x) 



1 



2nxX 



exp 



±(Y-Y')(y-y')-l-^- 
xX 8 Ci B 



{y + v'f 



y 



-i^-Bx{y + y') 



(6.3) 



where y=yi- y 2 , y' = y'i - y' 2 , Y = \{y t + y 2 ) Y' =A(v'i +J^) and x^-(x- x 1 ). 

If compared with the Green function derived in Sect.[y| Eq.(^^), the propagator Eq.( |6.3| ) contains an additional term 
proportional B, finite to the extent the distribution p(Q) is asymmetric. As we will see later, B creates the Lorentz 
force and plays in dynamics the role of an effective magnetic field. A finite Lorentz force in an Aharonov-Bohm array 
is not readily obvious: A classical Lore ntz force is absent since magnetic field is locally zero, whereas the quantum 
Aharonov-Bohm cross-section Eq. fl3.li] ) is left-right and <[> — > — $ symmetric and cannot explain B. Obviously, B 
and the associated force is directly related to the transverse momentum transfer considered in Sect. IB (see Eq.(A3)). 

In the limit of dense array with small typical flux, <f> — > 0, d AB — > 00, the effective field B reduces to to the 
macroscopic mean-field magnetic induction B — d AB < $ >. In this limit, the Aharonov-Bohm array model is 
equivalent the S— correlated field model Eq. (fLl|) with C~ 2 oc d AB < $ 2 >, in the external homogeneous magnetic field 
B. 

In general, however, ^-periodic B Eq.( |6.2| ) is very different from the mean-field expectations. Depending on the 
flux distribution function p($), the magnetic induction and B may be in any relation. For instance, the mean-field 
induction can be always compensated to zero by adding some properly oriented lines of flux = ^r&o, N = 1,2, 

In 



However, the added lines do not affect the effective magnetic field B seen by the particles (since sin ^27r^-J = 0). 
an Aharonov-Bohm array the Lorentz force may be finite even when the macroscopic magnetic induction is zero. 



A. Kinetic equation 



As before, in Sect.VB, consider the spatially uniform situation when the density matrix of the partially coherent 
wave at x = is of the form pi n (yi, ^2) = Po(y), y = y\ — yi- The density matrix p(y; x) of the wave at distance x 
can be found from Eqs.(4.4) and Eq.(|6.3[), 



p(y;x) = p (2/)exp 



.^yl-i^Bxy 



Again, as in Sect.VB. the density matrix obeys the following kinetic equation (compare with Eq.fc 



3 2tt 
V dx P + i ^^ VP + ^ ABP ^ ' Iai 



2C 



v 1 1 

—\v\ 

AB 



Performing the Fourier transform, one gets the equation for the distribution function n v in Eq.(5.!i): 

dn v eB dn v f 

v-^- + — ^ + I AB n v = . 

ox mc dip 



(6.4) 



The equation has the Boltzmann form and B enters kinetics as a magnetic field. The collision integral I AB has the 
form of Eq.(5.10) with the scattering rate 



W(4>) 



vX 1 



2tt£ 2 
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From here one can conclude that the relaxation is governed by an incoherent scattering by the flux lines: Indeed, W 
is proportional to the density of lines and the contribution of a line is given by the small limit of the Aharonov-Bohm 
cross-section Eq.( |3.11 ). Since even most dangerous small angle scatterin g fits this simple picture, it seems plausible 



that E q.(p.4[) can be generalized to arbitrary scattering angle using Eq.(3.11) as the probability scattering. Similar 
to Eq.(|5.14^ the kinetic equation reads 



where the collision integral 



g B din 

v ■ Vn„ tt^- + I A B n v> = , < w < 2ir 

mc dip 



2tt 

t - JL ( _ ^ 

1ab ^~ 2r AB J 27rsin 2 | ^ ^ ' 



with 



i n , / $ . 

— = -^(1-003 2.-) . (6.5) 



As before, the collision integral is a regular linear operator regardless its singular scattering-out term. Its action is 
defined by the following relation 

f im^ 1^1 im^ 

1 AB e ~ e 

AB 

From here, one sees that t ab has the meaning of the transport scattering time. 

Applicability of the Boltzmann equation requires the mean free path I ~ vt ab to be large on the scale of the wave 
length A. 

y ~A 2 d AB ^m 2 ^i^«l (6.6) 

If typically 3>/$o ~ 1, the density of the lines must not be too high: <i AB A 2 C 1. In case of lines with small $, the 
condition is milder: $ 2 e? AB A 2 <C 1. 

As an illustration, the Drude conductivity tensor can be readily derived from the Boltzmann equation Eq.( |S.5| ): 

a - = k* No * T " T+vh^ (6 - 7) 

a xy = (tt c T AB )a xx (6.8) 



where Nq is the density of states and fl c = — plays the role of the Larmor frequency. One sees that the Hall angle 
9 H , tan6» ff = <r xy /a xx , 



' sin 27r^- ' 

tanflff = -y^ "-t-y , (6.9) 

1 - cos2tt^ 

does not depend on the density of lines and parameters of the system. If the Aharonov-Bohm lines in array have same 
flux <I>, the Hall angle is simply 

t&nQ H = cot [ 7T— ] . (6.10) 



Counter intuition, the Hall effect is strongest at $ — > 0. 
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B. Landau quantization 



The Hall angle Eq.(3.10) in an Aharonov-Bohm array may be large which means that the particle orbit tends to 
be a circle, similar to Larmor orbits in an external magnetic field. The periodic classical motion is expected to be 
quantized (the Landau quantization) . To check this possibility, one should calculate the density of states of a particle 
moving in an Aharonov-Bohm array. In the quasiclassical region, when the energy of the particle E ^> Ml c , the 
problem can be generally solved by the method used in [ pT[ . In the present paper, only the most promising case of 
small flux array is considered. Not to repeat the calculation in [l4| , qualitative arguments which lead to the same 
result, are presented. 

In the linear in $ approximation, when scattering oc $ 2 can be neglected, the transverse momentum Eq.(3.12) is 
transfered to the particle, and it undergoes periodic motion with the angular frequency Q c = along the Larmor cir- 
cle with the radius Rl = v /VL C . In the small $-limit, the distinction between B and magnetic induction is immaterial, 
and the Bohr-Sommerfeld quantization condition can be formulated as the requirement of the flux quantization, 

$° rbit (E N ) = $ (N+l) , (6.11) 

where $ orblt (£") is the total flux encircled by an orbit of the particle with energy E. Due to the randomness in the 
flux line positions, the number of the encircled lines fluctuates, and so does the total flux. The total flux $ orblt = 
(<j>° rblt ) + ^$ is a random Gaussian variable which fluctuates around the average ($ orblt ) = BS{E), S(E) = -kR 2 l 
being the Larmor circle area. For the flux lines with uncorrelated positions, the total flux fluctuates with the variance 

((8<S>f)=d AB S(E)(&) • (6-12) 

Driven by the flux enclo sed by the orbit, the energy of the level fluctuates, En — (En) + SE. Neglecting the flux 
fluctuations in Eq.fl6.llp, $ orbit -> ($ orbit ), one gets the average energy of the N-th level (En) = hQ c (N + \)- To 



preserve the quantization Eq.(3.11), the level energy acquires a shift SE under the flux variations <5<j>: From the 
condition 

6$ + 6E n —\ En = , 
one gets the energy shift caused by the change of the flux (5$: 

se n = -m c — . (6.13) 







Combining Eq.(|6.13j) and Eq.Q6.12|), we see that the Landau levels acquire the Gaussian distribution form with variance 



(the width of the level) 



(6E N ) 2 ) = -— E (6.14) 



Physics here is similar to the inhomogeneous broadening: The Larmor circle sees different realizations of the random 
relief in different places of the "sample" , and the energy levels adjust their positions to the local conditions. 

A s alre ady noticed in |l4[ , it is rather unusual that importance of disorder (measured by the Landau level broadening 
Eq.( |6.14 )) increases with the energy E. The physical reason for this is that the larger is proportional to E the area 



under the Larmor circle, the bigger are the absolute fluctuations of the number of the Aharonov-Bohm lines encircled 
by the orbit and the flux fluctuation 

The densit y of s tates is given by the sum over the discrete levels. The position of the level is the Gaussian of the 
width in Eq.( |6.14 ) and centered at En = hfl c (N + |). The overlapping Landau levels create oscillating density of 



states. Similar to Ref. J14| , one applies the Poisson summation formula and finds the first harmonics of the oscillations: 
where the damping of the oscillations is controlled by 

2vr E 

7 
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Therefore, the Shubnikov-De-Haas oscillations in the Aharonov-Bohm array are strongly suppressed: even when 
O c r AB = tan Oh ^ 1 and the Larmor circling is well pronounced, quasiclassical quantization at high Landau levels 
E ^> fiVl c may not be seen because of the large damping 7. 
The damping parameter 7 can be also presented as 



1 



($ 2 



2 X 2 d^ (<f) 2 



> 



1 



2AV 



(6.15) 



where X (= fi/ \f2mE) is the wave length correspon ding to the energy E; the equality sign realizes in the case when 
the lines have equal flux. One concludes from Eq.(6.15) that as a necessary condition, the Larmor motion may be 
quantized only if the wave length exceed the distance between Aharonov-Bohm lines. As one could ex pect , this 
requirement tends to be complementary to the condition of applicability of the Boltzmann equation in Eq.(6.6). 



VII. CONCLUSIONS 



The main goal of the paper has been to understand specific features of the transport properties of a quantum charge 
subject to a random magnetic field, namely, the features related to the long range correlations of the gauge potential 
and the anomalous forward scattering. The non-perturbative method used in the paper to handle the divergence is 
based on the paraxial approximation to the stationary Schrodinger equation (Section nj). 

The paraxial theory of magnetic scattering is presented in Section |II A| . To show usage of the theory, it is applied 
to the Aharonov-Bohm line problem in Section HI. Being calculationally simple, the paraxial approximation proves 
to be rather efficient. The paraxial solution reproduces the small ang le a symptotics of the Aharonov-Bohm exact 
solution for the plane incident wave. The wave packet solution (see Eq.(3.9)) allows one to resolve an old controversy 
discussed in the Introduction concerning the transverse force exerted by the Aharonov-Bohm line. One sees that the 
angular distribution in the outgoing wave is indeed left-right asymmetric, so that there is a finite momentum transfer 
in the transverse direction. However, the asymmetry is concentrated within the angular width of the incident wave 
and it cannot be described in terms of the differential cr oss-se ction. For an arbitrary incident wave, the transverse 
momentum transfered to the charge can be found by Eq.(3.12). By comparison with the exact solution, the validity 
of this formula has been recently confirmed by Berry |}8| . 

The main result of the paper is the paraxial two-particle Green's function averaged with respect to the random 
magnetic field. In the paraxial theory, a stationary 2D problem becomes equivalent to a non-stationary ID one, and 
one can use the standard Feynman representation for the propagators. It turns out that the corresponding path 
integral can be evaluated exactly. The expressions for the Green's function is given by Eqs.([T^) and Eq.( |6.3| ) for the 
Gaussian random magnetic field and the Aharonov-Bohm array model, respectively. 

The paraxial two-particle Green's function solves the quantum problem of the near forward multiple scattering by 
random gauge potential: for a given incident wave, one is able to find correlators (rp(ri)ip*(r2)}, where averaging is 
performed with respect to the random field. To draw physical conclusions, two cases are analyzed: (i) (de)focusing of 
a coherent converging wave in the random magnetic field environment, and (ii) propagation of spatially homogeneous 
partially coherent beam. 

The non-local character of interaction with magnetic field is clearly seen from the analysis in Section |V A| of 
defocusing of a converging beam. The lost of coherence measured by defocusing is controlled by the size of the entire 
area "occupied " by the system that is the region where the wave function is finite: Indeed, it follows from Eq.(5.4) 
that the wave cannot be focussed if the random flux t hrea ding the area (aperture size) x (focus length) is of order of the 
flux quantum <!>o- Same conclusion follows from Eq.(5.7): having traveled a distance x, a perfect plane wave becomes 
a mixture of waves with the transverse momenta Ap y ~ hx/C 2 where the spatial coherence survives only within the 
region Ay ~ C 2 /x. Again, supporting the qualitative arguments presented in the Introduction, the coherence exists 
only within the area Aa; x Ay ~ C 2 ; the area is defined by the condition that the random flux is typically not bigger 
than the flux quantum $0. 

On the other hand, the ev oluti on i n the momentum space is rather ordinary: It is described by the usual Boltzmann 
equation derived in Sections VB and VIA. The only uncommon feature of the Boltzmann equation is that the collision 
integral being a well-defined operator nevertheless cannot be split in the scattering -in and -out regular parts: An 
attempt of the split produces two singular pieces, the two infinities canceling when combined. In the diagrammatic 
language, this can be rephrased as the cancellation of the divergences in the self energy and the vertex correction as 
observed in Ref. Q. 

The condition of applicability of the paraxial approximation can be derived from Eq.(2.9): At the distance x, the 
typical transverse size w is of order min[wo, C 2 /x] where wq i s th e characteristic length 
For large enough x, w ~ C 2 /x. Substituting this value into Eq^^), one gets 



1 in the initial distribution. 
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One sees that the theory is applicable in the non-perturbative region x 3> C in the paraxial limit X C. 

In the paraxial picture, the particle always moves in (almost) same direction. Therefore, any effect related to the 
Anderson localization is beyond the paraxial approximation. Although the Boltzmann equation allows for large angle 
scattering events, it is, of course, unable to describe the quantum localization either. Localization in random magnetic 
field remains a controversial issue, see Refs. Q and Jl9| . 

Two models of the random field has been considered in the paper: the Gaussian random magnetic field with zero 
average and the array of Aha rono v-Boh m m agnetic flux lines with arbitrary distribution of the line fluxes. Comparing 



the Green's function in Eqs. (5.2), and ( |6.3| ) (with B — ► 0), one sees that the models are paraxially equivalent. (For 



the Gaussian model, Eq.fl5.2|) has been derived assuming that the field is zero on average, (b) — 0. In a more general 



case, the Gaussian model with a finite magnetic induction B = (b) has the Green's function of the form in Eq.(3.S) 
with B substituted for B). 

The origin of the effective magnetic field in the Aharonov-Bohm array can be traced back to the left-right asymmetry 



in t he sc attering by an isolated Aharonov-Bohm line (see Section [II). One sees that the transverse momenta Ap± 
Eq.(|3~l~2|) gained as a result of collisions with Aharonov-Bohm lines, add together giving rise to the Lorentz force and 



the effective magnetic field B Eq.(3.2). Since any integer flux can be gauged out, Ap± and B are periodic functions of 



the fluxes. Note that the magnetic induction B, which by definition equals to ($), and the effective field B Eq.(5.2) 



are same quantities only if lines flux is small. Generally, they may be in any relation. In particular, one may have 

4 ' 2 



finite B ^ even if the magnetic induction B = 0: A complex of 3 lines with the fluxes: (+■%■, — %0 gives an 



example. 



In Section VI A , the kinetic equation for a charge in an Aharonov-Bohm line array is solved to find the Drude 



conductivity tensor Eq. (6. 7-3.5) and the Hall angle Eq.(p.9|). These results for the array of lines are in agreement 
with Ref. |Q where the Hall effect due to a single line has been studied. As discussed in |pH|| , the Aharonov-Bohm 
periodicity, i?$ + $ ~ i?$ combined with the time reversal symmetry, Bq, = — £>-$, requires that the half- integer flux 
lines, <f> = y do not generate any Lorentz force. The Abrikosov vortex carries the flux \ ■ — , therefore, does not 
exert the transverse force (of course, only in the limit when the particle wave length much exceeds the vortex size). 
As noticed before, this is in a qualitative agreement with the experimental observation of the reduced Hall effect 
exhibited by 2D electrons in the magnetic field of the Abrikosov vortices [[7p|JfO| . 

In an array of Aharonov-Bohm lines with small fluxes $ <C $o, the transverse momentum transfer oc $ (the 
"Lorentz" force exerted by B) is more efficient than relaxation of momentum the rate of which is oc $ 2 . Therefore, 

• ; ' TAR ' 



the "Lorentz" force significantly bends the p article trajectory leading to a large Hall angle Eq.(6.I0) and, probably. 



to Landau quantization. As shown in Section VI B , these expectation are almost never met. Even if the Hall angle is 
large and the Larmor circling is well pronounced, the Landau levels are very broad. The inhomogeneous in its nature 
broadening is due to the fluctuations in the total flux threading the Larmor orbit. 

Summarizing, the analysis of the forward scattering anomalous scattering in a random magnetic field (Gaussian 
and the Aharonov-Bohm array) has been presented using the paraxial approximation to the stationary Schrodinger 
equation. The gauge-invariant two-particle Green's function has been found by exact evaluation of the corresponding 
Feynman integral. The propagation of coherent (defocusing) and incoherent (Boltzmann equation) waves has been 
analyzed, as well as Landau quantization in the Aharonov-Bohm array. 
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APPENDIX A: THE TRANSVERSE FORCE 

The transverse force that is the transfered momentum in the direction perpendicular to the velocity, can be easily 
found in the paraxial theory. The expectation value of the transverse momentum in the outgoing wave is 

oo 

(p y )out = - J dyiP* out (x,y)—tp out (x,y) (Al) 

— OO 



23 



To calculate this integral, one exploits the fact that the free propagation at x > conserves the momentum, and, 
therefore, the average in Eq.(Al) does not depend on x > 0. The goal now is to present it as an integral at x = 
where the i/wt in Eq^J^l) is simplest possible; it cannot be done directly because of the eikonal discontinuity in i/wt 
at y = 

First, note that t/j on t{x,y) in Eq.(2.19), is a well-behaved function of y at any finite x > 0. Thus, the derivate 
in Eq.(Al) can be safely taken as the limit: -§^f{y) — ^- {f(y + if) — f(y — rj)) , r\ — > 0, and Eq. QAl| ) transforms to 



(Py)out 



4// 



(Pout (v)~ Pout (-»?)) , rj^Q 



(A2) 



where 



Poutiv) 



dyip? n (x,y-)ip in (x,y + ) , y± = y ± -77 



Substituting everywhere the subscript out — * m, one gets the transverse momentum in the incoming wave (p y )i n 
expressed via the corresponding Pi n (ii}- 

Identically, P ou t{rj) — (e^ v P y ). Since P out is the expectation value of the conserving variable e %v ^ y , its value does 
not depend on x. Choose the point x = +0 where Vwt is given by Eq.(|]l]) to evaluate P ou t- After simple calculations, 



PouM = PM + H(e 2 ^ -1) h/^ 



where 



Now, we see that the limit in Eq.(A2) is well defined, and we get for the transfered momentum Ap y = (p. 



y) out 



(P. 



y/in 



Ap y = ?i|i/> in (0)| 2 sin27r-J- 



(A3) 



This the final result for the transfered momentum. 
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From Eqs.( 2.17) combined with ( [2, 18] ) and (2.22), one gets Eq.(2.2J ) with the r.h.s. multiplied by exp[- 

This immaterial global phase factor is omitted in Eq.( ^2,23| ). 
. As discussed earlier in Sect. [I A, the overall phase of ipout does not have physicall meaning. Without changing any physics, 
ipout is multiplied by by exp[i7r<E>] for convenie nce . 
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Differentiating the l.h.s. of Eq.(4.6) with respect to x, one observes that the derivative vanishes by virtue of Eq.(fc 
means that the l.h.s. is x independent at x > x' and can be calculated in the limit x — > x' . In this limit G R (r, r')6{ 
5{y — y') and Eq.(|4.6j) follows immediately. 

Usually, when the scattering time r is finite, the initial distribution no decays as no exp[— t/r]. E.g. in the simplest case of 
the isotropic scattering, distribution function has the structure n(x) = no exp[— t/r] + (1 — cxp[— t/r] < no >). Here, the 
distribution smoothly transforms to the initial at t — > 0, and for t <C r scattering can be described by the perturbation 
theory. 

One sees from Eq.(5.8) that the eigenfunctions, \yo), I\yo) = Vo\yo) are 5-functions in the y 
S(y — yp ), and therefore, in the (^-representation, (<p\yo) = e r ~*~ . Then, /e 1- * - = (v/2C 



'M = 



Eq.(5.11) 



representation, 
yo e , and one comes to 
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